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Abstract. The behaviour of the toroidal and meridional components of the solar large-scale magnetic field 
and linear Alfven and Rossby waves during solar activity cycles and bi-annual time periods are theoretically 
investigated in this work. The shear of the toroidal and meridional components of velocity is taken into account. 
We consider the case of periodical velocity shear with bi-annual oscillation period, hereinafter the Velocity Shear 
Quasi Bi-annual Oscillations (VSQBO). The large-scale magnetic field toroidal and meridional components are 
obtained as harmonic functions of the time and the oscillation phase difference between them is equal to ir/2. The 
sign reversal of the magnetic field toroidal and meridional components is studied. The numerical simulations 
show that in the case of significant VSQBO amplitude values the toroidal or meridional component of the 
large-scale magnetic field reverses its sign three times in one of the hemispheres (northern or southern) of the 
Sun, during the solar activity cycle 23. According to our results the appearance of velocity shear oscillations 
leads to the modulation of the magnetic field 22-year period oscillations by the bi-annual ones. The presented 
model is applicable for investigation of the magnetic field evolution at the base of convection zone as well as for 
understanding the magnetic field properties in the upper solar atmosphere (the magnetic canopy of the solar 
chromosphere Roberts (1996)). 

The excitement of the linear Rossby and Alfven waves in the shear layer at the base of the convection zone 
is also considered. The periodical impulsive growth is characteristic to the energy density of Alfven and Rossby 
waves and they propagate as localized in time powerful pulses. Properties of these pulses are governed by the 
direction of waves phase velocity and the initial phase of VSQBO, and such behavior of waves well explains 
mechanisms of the solar flare excitement and activity. 
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1. Introduction 



The large-scale magnetic field longitudinal and latitudinal distribution and sign reversal during 
the solar cycle 23 have been investigated by many authors (de Toma et al., 2000; Benevolenskaya 
et al., 1999; Benevolenskaya, 1996; Benevolenskaya and Makarov, 1992; Bai, 1995; Bumba, 1989; 
Gaizauskas et al., 1983; Vitinskij, 1982; Ambroz, 1973). The recent observations of the p-mode 
frequency splitting and results of corresponding helioseimic inversions of the solar internal 
rotation show the appearance of two different layers with strongly marked shear flow (e.g. 
Kosovichev and Schou, 1997; Sekii, 1998). One of these layers is located at the base of convection 
zone and the differential rotation in this region, as is well known, plays an important role in 
the solar dynamo mechanism. The second one is located just below the photosphere. It is also 
reasonable to think that layer effects on the distribution of the global magnetic field at the 
solar surface and takes place in the formation of canopy of magnetic field in the upper solar 
atmosphere Roberts (1996). In other words the differential rotation takes an imposing part in 
the formation and evolution of the large-scale solar magnetic field. 

In this work we present the ideal magnetohydrodynamic (MHD) model explaining three-time 
sign reversal of the large-scale magnetic field at the beginning of cycle 23. It is known from the 
theory of the solar dynamo Parker (1979), that the toroidal and meridional components of the 
large-scale magnetic field oscillate with constant phase difference of n/2: the toroidal magnetic 
field reaches its maximum value when the meridional one is minimum or vice versa. In our 
model we obtain toroidal and meridional magnetic field components as harmonic functions 
of the time, when the differential rotation and its time oscillation are taken into account. 
Hereinafter in the work, the periodic time variation of the velocity shear with 2-year period is 
denoted as the Velocity Shear Quasi Bi-annual Oscillations (VSQBO). We present the result 
of numerical simulations showing that our model can explain one-time and three-time sign 
reversals of the large-scale magnetic field components during the activity cycle 23 and, in 
general, the modulation of activity cycle 23 by oscillations with a period of about 2 years. 
The behaviour of the large-scale magnetic field strongly depends on the amplitude and initial 
phase of VSQBO. The similar model was developed by Pataraya et al. (1997) using VSQBO. 
In that case the toroidal component of the magnetic field was a linear function of time and the 
meridional one was a constant. The structure of the large-scale magnetic field cyclic oscillations 
is also considered in the several works (Benevolenskaya, 1994, 1995, 1998; Akioka et al., 1987). 

We analytically investigate the possible mechanisms of the formation and evolution of the 
active regions and the corresponding excitement of the solar flares. In this sense, excitement 
of solar flares and their periodical activity are well known from several observations (e.g. see 
Bai, 1987; Rieger et al., 1984). In this work we investigate the propagation of the linear Rossby 
and Alfven waves in the environment with the steady and periodically oscillating velocity shear. 
Additionally, toroidal and meridional components of the magnetic field change sign one or three 
times during the activity cycle 23. We consider the time evolution of these waves. The numerical 
simulations show that the hydrodynamic and the magnetic energy densities oscillate with short 
periods and the waves propagate as localized in time pulses. These results are in agreement 
with the observed periodicity of solar flare activity and this allows one to think that our model 
is applicable for the explanation of solar flares excitement and activity. 



2. Basic Equations 



In this work we investigate the behaviour of the large-scale solar magnetic field in the layer, 
with the differential rotation and excitement of linear waves in this layer. As it is evident from 
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helioseismic observations there are two kinds of such layers in the solar interior. One of them 
is at the base of the convection zone and the other is just below the photosphere. These layers 
represent the areas bounded by two concentric spherical surfaces. The calculations are done in 
the co-rotating Cartesian frame, xoy plane is tangential to the bottom boundary of the layer. 
The x axes is directed eastward and it represents the toroidal direction, the y axes is directed 
northward and it represents the meridional direction and the z axes is directed outward to the 
upper solar atmosphere and it determines the radial direction. We use the ideal MHD equations 
for incompressible medium (Priest, 1982). So, for given latitude we have a set of differential 
equations in the following form: 

Dt dx y dy' K J 

(2) 

+ (fo + Py)V y , (3) 
-(fo + Py)V x , (4) 

(5) 

ux uy 

where, 

D d d d 

Dt = di + 'Hi + V (6) 

In these equations V x , B x are the toroidal and V y , B y are the meridional components of the 
velocity and the magnetic field, respectively; p = const and p are density and pressure; / = 
2Q sin 9 is Coriolis parameter; 

" = ^l (Sisin9) ' (7) 

where, Q is the angular rotation frequency of the sun; it!© is the Solar radius and 9 is the latitude 
corresponding to the origin of the local coordinate frame. The positive values of 9 correspond 
to the northern hemisphere and the negative ones denote latitudes in the southern hemisphere. 
The other symbols have their usual meanings. 

It is well known from the solar models that the rotation frequency of the Sun is a function 
of latitude £l = £1(6) and it could be well approximated as (Brown et al., 1989; Howard and 
Harvey, 1970): 

Q = Q (1- A lS m 2 6 - A 2 sin 4 6), (8) 

where, Qq, A\ and A 2 are constants. It should be noted that by using eq. (7) we obtain magnitude 
of the latitude, where (5 is equal to zero. After substituting A\ and A 2 values we have (3 — at 
the following latitudes: 

d[ N) w 74°33', 9[ S) w -74°33' (9) 

As we have mentioned above our aim is to study the excitement of the linear Rossby and 
Alfven waves in the layer where steady velocity shear and VSQBO appear. But before we turn 
to this matter let us determine the mathematical form of quantities satisfying the (l)-(5) set 
of equations in an unperturbed state of environment. 
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3. The Sign Reversal of the Magnetic Field Components 

We consider the differential rotation of the Sun with respect to the co-rotating Cartesian frame 
of reference described above. In this case we have the toroidal component of the unperturbed 
velocity in the following form: 

V x0 = yV xy , (10) 
where, V xy is the shear in the y direction of the toroidal flow and it has the form: 

V xy = V + V xn cos(uj n t + <p ) = Vf 2 , (11) 

The first term on the right side of eq. (11), V is the permanent shear of the toroidal flow in the 
considered layer and it conveys the steady differential rotation of the Sun: 

V = -^(n cos 9) = -Qq sin 9(1 + 2A 1 + (AA 2 - 3A ± ) sin 2 9 - 5A 2 sin 4 9). (12) 
d9 

This expression was presented in Pataraya and Zaqarashvili (1995) and Gigolashvili et al. 
(1995). [In these works the excitement of linear waves were investigated, when the toroidal 
component of the unperturbed magnetic field is a linear function of time.] 

By using eq. (12) we can obtain the latitude value, where the steady velocity shear of the 
toroidal flow V has a maximum by solving equation dV/d9 = 0, and after substituting of A 1 
and A 2 values we obtain: 

9f ] w 55°9', 9 { 2 S) « -55°9' (13) 

The other term on the right side of eq. (11) determines the velocity shear oscillation with 
an arbitrary period and in the case of 2 year oscillation period it demonstrates the VSQBO. In 
this term u n is the frequency of the velocity shear oscillations: 

= — , (14) 

nt 

where t ~ 3.16 • 10 7 sec= 1 year. nt is the oscillation period of the velocity shear. 
We have the meridional component of the unperturbed velocity in the form: 

V y0 = xV yx , (15) 

where V yx is the shear of the meridional flow in the x direction: 

o 2 

V yx = -^/ 2 , (16) 

According to eqs. (10)-(12), (15)-(16) and the equation of induction we obtain expressions 
for the toroidal B x0 and meridional B y0 components of the unperturbed magnetic field: 

B x0 = -B xm sin(fiit /i), (17) 

B y o = B ym cos(fiitoA), (18) 

where, Qi = 2n/22t is the frequency of low-frequency oscillations (with the period 22 years) 
of the unperturbed toroidal and meridional components of the solar magnetic field; B xm = 
B sin 9 = V ■ Bym/Qi here, B = const 

/i = t - 1 + a n [sin( — r + ip Q ) - sin(^ )], (19) 
n 
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where r = t/t and 



/ 2 = l + % 1 cos(-r + v.o), (20) 
V n 



lr " (21) 



u n t V 

As we can see from eqs. (17) and (18) the toroidal and meridional components of the magnetic 
field are harmonic functions of the time and they oscillate with a constant phase difference of 
7r/2. We carried out the numerical simulations for the case of n = 2 in eq. (14). In other words, 
we consider the VSQBO. We studied the time evolution of the magnetic field components 
defining the VSQBO amplitude in two different ways: 

a) Using in eqs. (11) and (21) the following expression of the V xn : 

V xn = V sm 2 9cos9, (22) 

one has 

Vq sin 9 cos 6 

a2 ~~7iQol + 2A 1 + (AA 2 -3A 1 )sm 2 9-5A 2 sm 4 9' ^ ' 

Analysing the results of numerical simulations we found that the low-frequency (with 22-year 
period) oscillation of the large-scale magnetic field components is modulated by the high- 
frequency (with bi-annual period) one. The modulation rate depends on the latitude 9 and 
the magnitude of a 2 . Particularly, the significant modulation of the magnetic field oscillation 
takes place, when n\a 2 \ > 1 and the three-time sign reversal of the magnetic field toroidal 
component occurs at the beginning of the activity cycle. Figure 1 shows the time evolution 
of the normalized toroidal B xQ /B xm (solid lines) and the meridional By /B ym (dashed lines) 
components. Figure la corresponds to the latitude value, where the steady velocity shear of the 
toroidal flow has maximum 9 = 9 2 N ^ in the northern hemisphere and Figure lb corresponds to 
the same latitude in the southern hemisphere 9 = 9^ (see eq. (13)). In this case we assume 
that the amplitude of the VSQBO V xn is defined by eq. (22) and its initial phase ipo = 7r. These 
results show that at the middle latitudes the amplitude of VSQBO is large enough to induce 
the significant modulation of the large-scale magnetic field cyclic oscillations by the bi-annual 
ones in both hemispheres of the Sun. This modulation leads to the magnetic field toroidal 
component three-time sign reversal within the first two years of the activity cycle 23 in the 
northern hemisphere (see the encircled area 1 in Figure la), while in the southern one only 
one-time sign reversal occurs during the same time interval. As one can see from the presented 
Figures, when tp = n at the solar activity minimum the three-time sign reversal of the magnetic 
field toroidal component appears only in the northern hemisphere. We have a similar picture 
in the southern hemisphere in the time interval between 11 and 13 years from the beginning of 
the cycle 23. The time behaviour of the magnetic field toroidal component is inversed within 
this time interval, i.e. now three-time sign reversal occurs only in the southern hemisphere 
(see the encircled area 1 in Figure lb). In other words, the modulation of the magnetic field 
components oscillation is not synchronized between different hemispheres and, because of this, 
we have the asymmerty in the three-time sign reversal occurrence. This kind of behaviour is 
characteristic of the magnetic field toroidal component for most latitudes except the areas very 
close to the equatorial plane and pole. We will address this in more detail later. With regard to 
the magnetic field meridional component (the dashed line curves), its low-frequency oscillation 
is also strongly modulated by the bi-annual one in both hemispheres of the Sun. In spite of 
this, the three-time sign reversal of the meridional field component disappears for most of the 
latitude values, although this effect can almost be achieved only close to the latitudes where the 
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steady velocity shear is maximum (eq. (13)) and this effect occurs in the northern hemisphere 
in the time interval 6 < r < 8 years (see the encircled area 2 in Figure la) and in the southern 
one in the time interval 5 < r < 7 years (see the encircled area 2 in Figure lb). It is important 
to note that, if we perform calculations for the case of ipo = 0, then we obtain an inversed 
picture, i.e., in this case the three-time sign reversal of the magnetic field toroidal component 
occurs only in the southern hemisphere, within the first two years of activity cycle. 

As we mentioned above the situation changes in the equatorial and polar areas. The time 
evolution of the magnetic field components in these areas is very similar in both hemispheres, 
we therefore present the numerical results regarding the northern hemisphere. In Figures lc and 
Id the curves for latitudes close to the equatorial plane (9 = 5°) and to the pole (9 = 85°) in 
the northern hemisphere are presented, respectively. The behaviour of the large-scale magnetic 
field in the equatorial and polar areas is different from that considered above. In these areas 
the amplitude of the VSQBO is small and, because of this, 7r|a2 1 becomes < 1. As one can see 
in Figures lc and Id, in this case the very weak modulation of the magnetic field low-frequency 
oscillation occurs and the above mentioned sign reversal does not appear. Finally, exactly on 
the equatorial plane (9 = 0°) or on the poles (9 = ±90°) the modulation of the cyclic oscillation 
vanishes completely. 

b) We have also performed the numerical simulations for the case when the amplitude of the 
VSQBO V xn is defined as: 

V xn = -U Vsm9(5(l - (Lu n t )- 1 ) + sm9), (24) 
and in this case we obtain a 2 as follows: 

U sin 9(5(1 -7T- 1 )+ sin 9) 
a 2 = , (25) 

71 

where U and 5 are constants. As in the first case, we present in Figure 2 the curves showing 
the time dependence of the normalized toroidal and meridional magnetic field components, 
when the VSQBO amplitude is defined by eq. (24), its initial phase is ipo — 7r, 5 — 1.4 and 
Uq = 1.8. Figure 2a corresponds to the latitude value in the northern hemisphere, 9 = 9^ N \ 

(S) 

and Figure 2b to that in the southern one, 9 = 9 2 ■ These results show, that unlike the 
previous case, when 5 = 1.4 the strongly marked modulation of the magnetic field component 
low-frequency oscillation appears (i.e. n\a 2 \ > 1) only in the northern hemisphere, while in 
the southern one 7r|a2| < 1 and the modulation is weak. Hence, the three-time sign reversal 
of the magnetic field toroidal component within first two years of the activity cycle occurs 
only in the northern hemisphere (see the encircled area 1 in Figure 2a), while in the southern 
hemisphere the modulation is very weak and thus the three-time sign reversal is absent in 
the whole hemisphere. Unlike the previous case, the asymmetry in the magnetic field toroidal 
component three-time sign reversal is not due to different phases of modulation of the low- 
frequency oscillations, but now the strong modulation occurs only in the northern hemisphere. 
As it is shown in Figure 2a the oscillation of the meridional component is modulated too, but in 
spite of this its three-time sign reversal does not occur at most of latitude values. It is important 
to note that, we can obtain an inversed picture in the case, when 5 = —1.4 and (p = 0. i.e. 
in this case we have the strong modulation of the magnetic field cyclic oscillation only in the 
southern hemisphere. 

Finally, the behaviour of the magnetic field components close to the equatorial plane is 
similar to that shown in the previous case. As one can see from Figure 2c, the above mentioned 
modulation is very weak and thus the sign reversal effect does not occur near the equatorial 
plane. Although, from eq. (25) one can see that | «2 1 increases with latitude 9 and when S > it 
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is maximum on the northern pole and when 5 < it is maximum on the southern pole. Hence, 
in these cases the three-time sign reversal the magnetic field toroidal component appears at 
one of the polar areas too (see the encircled area 1 in Figure 2d). Moreover, in this case the 
three-time sign reversal of the meridional component can also be achieved (see the encircled 
area 2 in Figure 2d). 



4. The Evolution of Rossby and Alfven Waves 

In the previous section we studied (l)-(5) a set of differential equations in the unperturbed 
state and we determined the expressions for the toroidal and meridional components of the 
unperturbed flow velocity and magnetic field. Let us now consider the excitement of the linear 
waves in the incompressible medium located in the solar internal shear layer. In the unperturbed 
state we take into account the shear and VSQBO (eqs. (10) and (15)). Additionally, the 
components of the magnetic field are oscillating with period of 22 years. At the same time 
this low-frequency oscillation is modulated by the high-frequency oscillation of the magnetic 
field with period equal to that of VSQBO. 

We investigated equations (l)-(5) by using the linear perturbation theory. So, all quantities 
are presented in the form ip = ipo+ip and we obtain the following set of the linearized differential 
equations : 

f^ = 0, (26) 
ox dy 

D o B 'x ( „ d „ d 



D t 



B ^ + B y o£)v; + V yx B' x , (28) 



A) ( dv' x dV\ 1 / d , d\(dB' x dB y 



Dot \ dy dx J Airp 

where, 



IH +,/ 4+ i '4 (30) 

As we can see from eqs. (10) -(15) both V x o and V y o depend on the time and spatial coordinates. 
On the other hand the magnetic field components are harmonic functions of time (eqs. 17-18). 
Because of this, in the equations we have temporal as well as the spatial inhomogeneity. In this 
case the modal analysis is not applicable for the investigation of the perturbed quantities. 

We studied eqs. (26)- (29) by using the non-modal perturbation theory which is considered 
by many authors (Chagelishvili et al., 1993; Lominadze et al., 1988; Farrell, 1987; Marcus 
and Press, 1977; Taylor, 1931). According to non- modal analysis, when we consider the linear 
velocity shear profile, it is possible to find out the coordinate transformation, which in the 
linearized equations, leads to the changing of spatial inhomogeneity by a temporal one. Using 
this coordinate transformation allows one to investigate the temporal evolution of the spatial 
Fourier harmonics and in this case all of the perturbed functions have the following form: 

V>' = i)\exp{iK x x + iK y y}, (31) 

where, 

K x = k xl F u (32) 
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K y = ^-k xl F 2 , (33) 

are the toroidal and meridional components of the wavevector, respectively. In expressions (32) 
and (33): 

F 1 = CK 21 sin^fo/i) + cos^^oA) (34) 
F 2 = - sin(fit /i) + CK 21 cos(m /i), (35) 

where, C = Q\/V and K 2i = k y \jk x \ In these equations k x \ and k y \ are constants. According 
to expressions (26)-(29) and (31)- (33) we obtain the set of equations for perturbed functions 
in the following form: 

K x B xl + K y B yl = 0, (36) 
K x V xl + KyV y i = 0, (37) 

(W " X = iy/Wpu A V yl - ^V yx B yl , (38) 



dt v r y K x 



where, 



Kl = Kl + tf* (40) 



y 

and ^ is the frequency of Alfven waves and it has the following form: 

= K^ + KyByO = j^ Ck (41) 

where, -E^™ is the amplitude of the magnetic field toroidal component (see eq. (17)). First, we 
consider the case when k yl — 0. In this case u A = and propagation of only Rossby waves is 

( Ft) 

possible and we obtain the analytical solution V y \ = V y \ of eq. (39) in the following form: 

v - = v ^ - + (c-^SW) M- ^-"- ^ ^[(--D^-^/Ol} 

(42) 

Since all perturbed quantities are assumed to be ~ exp(iip) and the oscillation frequency uj = 
—d(p/dt we obtain the dispersion relation of Rossby waves: 

= — — T S ^f "Alo : f u + ^1^1/2 sin(fiitoA) + yk x iVtif 2 cos(fii* /i), (43) 

fcxi(l + (C~ 2 - 1) sin (fiito/i)) 

and at last, by using of eqs. (37)and (42) we have the expression of Rossby waves hydrodynamic 
energy density W R = p(\V^\ 2 + \ Vyf > \ 2 )/2 in the following form: 

Wr = « ! (44) 

2 l + (C- 2 -l)sin 2 (!V /i) 

In Figure 3 the time dependence of the W R normalized by its initial value is presented. Figures 
3a and 3b correspond to the northern (6 = 6 2 N ^) and southern (6 = O^) hemispheres, 
respectively, when the amplitude of the VSQBO is defined by eq. (22) and its initial phase 
ipo = 71. As we can see from these figures, when in the medium propagate only Rossby waves 
their hydrodynamic energy periodically increases in short time period during the solar activity 
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cycle, which is demonstrated by peaks in the plots. The time moments of the wave energy 
density peak appearance coincide with those when the unperturbed magnetic field toroidal 
component changes its sign. As we have mentioned above, in this case occurs three-time sign 
reversal of the toroidal component and this fact reveals itself in the three peaks at the beginning 
of the activity cycle 23 (Figure 3a). The peaks denoted by "1" in Figures 3a and 3b represent 
the the moment of sign reversal due to the mean low-frequency oscillations of the magnetic 
field toroidal component and the occurrence of peaks denoted by "2" and "3" is caused by 
the three-time sign reversal effect. Hence, when one-time sign reversal occurs one has only 
the single peaks. In the moments of sign reversal of the unperturbed magnetic field toroidal 
component the Rossby waves propagate as pulses localized in time. [It should be noted that by 
comparison of Figure 3a and 3b it is clearly seen that the asymmetry between the northern and 
southern hemispheres due to the fact that modulation of the magnetic field cyclic oscillation in 
the different hemispheres is asynchronous.] 

We have also performed the numerical simulations in the case when the amplitude V xn of 
the VSQBO is defined by eq. (24), its initial phase (p — 7r, S — 1.4 and U = 1.8. Resulted 
curves presented in Figure 4 showing the time dependence of the Rossby wave hydrodynamic 
energy density normalized by its initial value. Figures 4a and 4b correspond to the northern 
(9 = 9 2 ) and southern (9 = 9 2 ) hemispheres, respectively. In these Figures we have the 
peaks of energy density in the time moments when sign reversal of the unperturbed magnetic 
field toroidal component occurs. As we know from the above discussion in this case we have 
three-time sign reversal only in the northern hemisphere and this leads to the appearance of 
two additional peaks denoted by "2" and "3" beside the mean peak denoted by "1" in Figure 
4a. Furthermore, we have only single peaks in the southern hemisphere and this is a result of 
the disappearance of the strong modulation of the magnetic field cyclic oscillations by bi-annual 
ones. This property of waves determines the asymmetry of solar flare activity intensity within 
the solar hemispheres and a detailed discussion on this matter is presented below. 

It is important to note that, in general, in the case of k y \ = 0, the meridional component 
of the Rossby wave wavevector is K y = at the time moments when the sign reversal of the 
unperturbed magnetic field toroidal component occurs. This means that, in time moments when 
the Rossby waves propagate exactly in the toroidal direction the wave energy density increases 
and it has the form of a sharp peak. 

Now let us consider the case when k yl ^ 0. In this case, in the medium, propagation of both 
Rossby and Alfven waves is possible. We performed a numerical calculation of the (36)-(39) 
set of differential equations. In Figure 5 the normalized hydrodynamic and magnetic energy 
densities are presented as functions of time corresponding to 9 = 9 2 N ^ (Figures 5a) and 9 = 9 2 S ^ 
(Figures 5b) values of the latitude. One can see from these curves that wave hydrodynamic 

W H = P -{\V xl \* + \V yl \>) (45) 

and magnetic 

W M = l^ + l^l 2 (46) 

energy densities periodically increase during the short time intervals while in the initial mo- 
mentum we assumed that the hydrodynamic and magnetic energy densities are commensurable 
Wh(0)/Wm(0) ~ 1- The linear Alfven waves propagate as powerful pulses. The period of energy 
density oscillation depends on the k xl and k y \ magnitude. In other words the hydrodynamic 
and magnetic wave energy densities oscillate with period which depends on the direction of 
the wave phase velocity in the r = time moment. Figure 5 shows curves for k x \ = 10 -9 
cm -1 and K 2 i = 2. The calculations are carried out for values of environment parameters 
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satisfying those obtained from the standard solar models for the base of the solar convection 
zone. Particularly, the amplitude value of the unperturbed magnetic field toroidal component 
(eq. (17)) is determined assuming B = 10 5 G. Figure 6 shows the time variation of the 
hydrodynamic and the magnetic energy densities for the same values of parameters, but in this 
case the amplitude of the VSQBO is taken by eq. (24). As it is clear from the curves in Figure 5 
and Figure 6 the character of the wave temporal evolution is similar in both cases. This model 
well explains the mechanism of the solar flares excitement and activity. For the different values 
of k x i and k y \ the wave energy density oscillation with periods well describing the observed 
51-day (e.g. see Bai, 1987) as well as 154-day (e.g. see Rieger et al., 1984) periodicity in the 
occurrence of the solar flare activity may be obtained. Observational data (Bai, 1987) also 
show asymmetry in intensity of the flare activity between northern and southern hemispheres 
of the Sun. Comparing curves demonstrated in Figures 5a and 5b (Figures 6a and 6b) one 
can see that the amplitude of energy density oscillation in northern hemisphere is greater 
than that in southern hemisphere and the more intensive excitement of the solar flares in the 
northern hemisphere is expected. According to our results the flare activity is more intensive 
in the hemisphere where the thee-time sign reversal of the unperturbed magnetic field toroidal 
component occurs within the first two years of the activity cycle. 

As we have shown in section 2 in the unperturbed state the modulation of large-scale 
magnetic field low-frequency oscillation vanishes close to the polar and equatorial areas, when 
the VSQBO amplitude is determined by eq. (22). This leads to the decrease of the wave energy 
oscillation amplitude in mentioned areas. On the other hand, when we consider the case of the 
VSQBO amplitude by eq. (24) then in the polar area of the northern hemisphere the amplitude 
of the wave energy density oscillations remains to be significant. 

So, we have presented the model of solar flare excitement at the base of the solar convection 
zone as the reasonable mechanism of the solar flare activity origin. As we mentioned in the 
introduction our model is also applicable for the shear layer just below the photosphere. It is 
reasonable to think that the velocity shear effects on the large-scale magnetic field distribution 
at the solar surface and the excitement and the time evolution of Rossby and Alfven waves take 
part in the formation of the magnetic field structure in the upper solar atmosphere and in the 
heating of the solar corona. 

5. Conclusions 

In this work we investigated the behaviour of the large-scale magnetic field in the solar internal 
shear layers. The VSQBO is taken into account. By use of the numerical simulations we studied 
the behaviour of the magnetic field toroidal and meridional components. The excitement and 
time evolution of the linear Rossby and Alfven waves at the base of convection zone were also 
studied, when the VSQBO and modulation of the magnetic field components were taken into 
consideration. We obtained results for our analytical solutions and for the above mentioned 
numerical simulations, which can be formulated as follows: 

1) We obtained the expressions of the magnetic field toroidal and meridional components in 
the form of harmonical functions of time. 

2) We obtained the curves representing the time dependence of the magnetic field toroidal 
and meridional components in two cases: 

a) When the VSQBO amplitude is defined by eq. (22). In this shown in Figures la 

and lb, we have the strongly marked modulation of the large-scale magnetic field low-frequency 
(with period 22-years) oscillation by the high-frequency (with the 2-year period) ones at most 
latitudes of both northern and southern hemispheres of the Sun. When the initial phase of 



The Magnetic Field Sign Reversal and Evolution of Rossby and Alfven Waves in the Solar Interior 



11 



VSQBO (f = 7r, the modulation of the magnetic field components cyclic oscillation leads to 
three-time sign reversal of the toroidal component in the northern hemisphere, within the first 
two years after activity minimum. During the same time interval we have only one-time sign 
reversal of the magnetic field toroidal component in the southern hemisphere. This asymmetry, 
in occurrence of the three-time sign reversal effect between different hemispheres, is caused 
by asynchronous modulation phase in these hemispheres. The three-time sign reversal of the 
magnetic field meridional component does not happen at most latitude values. Exceptions are 
areas where the steady shear of the velocity toroidal component is maximal (see eq. (13)). In 
this case the three-time sign reversal of the meridional component may happen in time interval 
between 6 and 8 years after activity minimum in the northern hemisphere and between 5 and 
7 years in the southern hemisphere. 

The above mentioned effect of the modulation decreases in equatorial and polar areas (Fig- 
ures lc and Id). Because of this, in these areas three-time sign reversal of any magnetic field 
component becomes impossible. 

b) We also considered another case, when the VSQBO amplitude is defined by eq. (24). 
Corresponding curves are presented in Figure 2. In this case, when the initial phase of VSQBO 
ipo = 7r the strongly marked modulation of the magnetic field components low-frequency oscilla- 
tion appears only in the northern hemisphere (see Figure 2a), while in the southern hemisphere 
this modulation is very weak (Figure 2b). Because of this the three-time sign reversal of the 
field toroidal component occurs only in the northern hemisphere, within the first two years after 
activity minimum. This asymmetry is the result of the modulation range difference between 
solar hemispheres. 

Similar to the previous case the modulation range decreases in equatorial area, where the 
three-time sign reversal effect vanishes. But in this case the magnetic field cyclic oscillation 
modulation range is maximal in the northern polar area. Hence, in this area modulation remains 
significant. Moreover, near the northern pole the three-time sign reversal of the magnetic field 
meridional component becomes possible. 

Further, we considered the excitement of the linear Rossby and Alfven waves in the base 
of convection zone, when in unperturbed state VSQBO and corresponding modulation of the 
magnetic field cyclic oscillations are taken into account. We carried out the analytical and 
numerical investigation of the excitement and evolution in time of these waves and we obtained 
the following results: 

3) When in eqs. (34)-(35) k yl = 0, i.e. K 2 \ = in the environment propagate only Rossby 
waves. In this case it is possible to analytically solve the set of differential equations (36)-(39) 
and we have obtained a solution of the meridional component of the perturbed velocity (see eq. 
(42)). By using this solution we determined the expression of the Rossby wave hydrodynamical 
energy density temporal variation. 

4) We obtained the curves of the Rossby wave energy density time dependence. Curves 
in Figure 3 correspond to the definition of the VSQBO amplitude by eq. (22) and Figure 4 
corresponds to that defined by eq. (24). These figures show that the Rossby wave energy density 
increases during very short time intervals and such behaviour of the wave energy density are 
shown as sharp peaks. These peaks appear in the time moments when the sign reversal of 
the unperturbed magnetic field toroidal component occurs. Thus, in time moments when the 
three-time sign reversal of the magnetic field toroidal component happens we have three peaks 
instead of one. In the case of k y i = these time moments coincide with those of the wavevector 
meridional component K y — 0. This means that in these moments Rossby waves propagate 
exactly in the toroidal direction. 

5) We performed the numerical investigation of the set of differential equations (36)-(39) in 
case when k yi ^ 0. In this case propagation of both Rossby and Alfven waves is possible. We 
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obtained the curves of the wave hydrodynamic and magnetic energy density temporal variation. 
Figure 5 corresponds to the case "a" of the VSQBO amplitude and Figure 6 corresponds to 
case "b". These results clearly show oscillation of the wave hydrodynamic and magnetic energy 
densities. By choosing the k x \ and k y \ values one can obtain oscillation of wave hydrodynamical 
and magnetic energy densities with periods comparable with those obtained from observations 
of the solar flare activity (e.g see Bai, 1987; Rieger et al., 1984). These reasons allow us to 
think that the model we have presented here well describes the solar flare excitement and 
their periodical activity. Our results are in agreement with the observations of solar flare 
activity which show the difference of the flare activity intensity between northern and southern 
hemispheres of the Sun (Bai, 1987). 

The model presented here is applicable for investigation of the magnetic field variations and 
the linear wave evolution in the layers where the arbitrary period oscillations of the velocity 
shear, characteristic to Solar atmosphere, occur. In this sense it is important to study the role 
of the shear layer just below the photosphere in the distribution of the magnetic field at the 
solar surface. This model is also useful for the investigation of the excitement and evolution 
of waves close to the solar surface and their influence on the solar corona. Finally, taking into 
account the velocity shear oscillation with 11-year period, these work results could be used for 
investigation of solar p-mode properties during activity cycles. These very important tasks will 
be subjects of our forthcoming works. 
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Toroidal and Meridional Components 
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Figure 1. The normalized toroidal B x0 / B (m (solid line) and meridional B /B (dashed line) 
components of the solar large-scale magnetic field are plotted against x , time (Years). The 
VSQBO amplitude is defined by eq. (22) and its initial phase is (p =ji . a) the latitude = { 2 N) ; b) 

9 = 6 ( 2 S) ; c) the equatorial area 9 = 5° ; d) the polar area 9 = 85° . 
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Figure 2. The normalized toroidal B x0 /B xm (solid line) and meridional B /B (dashed line) 
components of the solar large-scale magnetic field are plotted against x , time (Years). The VSQBO 
amplitude is defined by eq. (24) and its initial phase is (p Q = K . a) the latitude 9 = ( 2 N) ; b) 6 = ; 
c) the equatorial area = 5° ; d) the polar area = 85° . 
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Figure 3. The curve of the normalized Rossby wave hydrodynamical energy density 
W R (t)/W R (0) as a function of t , time (Years). The VSQBO amplitude is defined by eq. (22) and 



its initial phase is <p = % . a) the latitude = 0^ 5 b) 0=0 



)(S) 



iiie nyuiuuynamiu cneigy i^eiisny ui rvussuy waves 



w R (o) 



4 - 



North 



2 - 



a 



T 



T 



8 12 

T.TIME (Years) 



16 



20 



A l\ , 



1 

24 



™*Q R xlO 3 



4 - 



2 - 



South 



-i t 1 r 

12 16 
X, TIME (Years) 



20 



~ I 
24 



Figure 4. The curve of the normalized Rossby wave hydrodynamical energy density 
W R (z)/W R (0) as a function of t , time (Years). The VSQBO amplitude is defined by eq. (24) and 

its initial phase is (p = K . a) the latitude 9 = ( 2 N) ; b) 9 = ( 2 S) . 
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of Linear Alfven Waves 





Figure 5. The normalized Alfven wave hydrodynamical W H (t)/W H (0) (solid line) and the 
magnetic W M (t)/W m (0) (dashed line) energy densities are plotted against T , time (Years), when 
in the initial momentum W H (0)/W M (0) ~ 1 . The VSQBO amplitude is defined by eq. (22) and its 
initial phase is (p = K . a) the latitude Q = ; b) 9 = 9 ( 2 S) . 
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Figure 6. The normalized Alfven wave hydrodynamical W H (t)/W H (0) (solid line) and the 
magnetic W M (t)/W M (0) (dashed line) energy densities are plotted against T, time (Years), when 
in the initial momentum W H (0)/W M ~l. The VSQBO amplitude is defined by eq. (24) and its 
initial phase is (p = K . a) the latitude = 6 ( 2 N) ; b) = ( 2 S) . 



